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ABSTRACT 

This thesis is a continuation of a series of studies on automatic 
throttle compensation systems begun in 1963 by Lieutenant Commander 
С. В. Bell. [1] Previous investigations have been concerned with flight 
in a wings level attitude. This study will consider the effects of tuming 
flight on the system. The system under consideration here is the Special- 
ties, Incorporated, Automatic Power Compensator, (APC), an automatic 
throttle in use in the F-8 and F-4 aircraft. It utilizes feedbacks of angle 
of attack and normal acceleration. И was found important to differentiate 
between spurious perturbations and those perturbations from steady state 
introduced due to the aircraft entering a turn. 

This study was conducted by Lieutenant Richard D. Hartman at the 
United States Naval Postgraduate School, Monterey, California. The 
assistance given by Professor E. J. Andrews: is greatly appreciated by 


the writer. 
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L Introduction. 

This study is one of a series being conducted at the United States 
Naval Postgraduate School to evaluate the use of automatic throttle com- 
pensation systems in modern high performance jet aircraft. The need for 
automatic throttle compensation was discussed thoroughly in the previous 
investigations. [1 " 2 | Briefly, the slope of the thrust required versus air- 
speed curve is negative at the approach airspeed for modern jet aircraft. 
This means that a loss of airspeed may rapidly develop into an uncontrolled 
sink rate. Aircraft with this characteristic are said to have negative static 
airspeed stability. Positive airspeed stability must be provided by proper 
use of the throttle. An automatic throttle compensation system relieves 
the pilot of the task of providing this airspeed stability. With an automatic 
throttle system keeping a constant airspeed, the pilot can concentrate on 
maintaining his altitude or glide slope by use of the elevator. 

There are various combinations of errors in airspeed, pitch attitude, 
angle of attack, and normal acceleration that can be used as inputs to an 
automatic throttle device. Direct use of these feedbacks and their deriva- 
tives and integrals provide the final throttle command to the engine for the 
desired change in thrust. The system under consideration here utilizes 
feedback of angle of attack and normal acceleration. A block diagram of 
this system is shown in Figure 1. The transfer functions, also given in 
Figure 1, are those from Reference 3. 

д Discussion. 


Previous studies conducted at the U. S. Naval Postgraduate School 


considered flight in a wings level attitude. It is the intent of this study to 
extend the analysis to include turning flight and to investigate the possible 
additional inputs to the basic automatic throttle system that may be required. 
The aircraft parameters used are for the F-8, an operational U.S. Navy swept- 
wing fighter. Previously, the equations involved were simulated on an 

analog computer and a obtained in this. manner. For this investiga- 
tion, use was made of the 1604 Control Data Corporation digital computer 
atthe U. 5. Naval Postgraduate School. A program titled "Runge-Kutta 
Solution of Ordinary Differential Equations with Built-in Input and Output 
Routines" was used to provide solutions for the equations. 

Since the digital computer had not been used in References 1 and 2, 
the response of the basic airframe in a wings level attitude was first exam- 
ined. This preliminary investigation was conducted with and without the 
automatic throttle compensator in use. When the results of this investiga- 
tion showed close agreement with other results, turning flight was investi- 
gated. 

This analysis is confined to the longitudinal mode of operation and 
the linearized three-degree-of-freedom, body axes, equations of motion 
derived from small perturbation theory are utilized. These equations as 
derived in Reference 4 are: 
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The values for the stability derivatives in the above equations were 
obtained from Reference 3 and are given in Table 1. These values are for 
the F-8 aircraft in landing configuration at approach speed. Several of these 
parameters are functions of Cr, and C p and would change values if found 
advisable to change the aircraft angle of attack while in a turn. Utilizing 
the values of the parameters and their basic definitions, they are developed 
in Appendix I as functions of Cry, and C p. 

Curves of C; and Cp versus angle of attack in Reference 5 were used 
to obtain the plot of Cy, versus Cp shown in Figure 2. An equation of this 
plot is derived in Appendix II and found to be 

CD = .955 - 2.12(С1) + 1.44(С1)2 

In order to write a program for the digital computer a flow graph of 
the system to be analyzed is first developed. Such a flow graph for the air- 
frame and automatic throttle combined is shown in Figure 3. The program 
for the digital computer developed from this flow graph is shown in Figure 4. 
Minor variations in this program and the use of initial conditions allows for 
an analysis of the system under a great variety of conditions. 


3. Flight with Wings Level. 
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Aircraft Without Automatic Throttle 

By eliminating the change of thrust inputs to the moment and X-force 
equations, the automatic throttle will be disconnected from the aircraft 
equations of motion. This was easily effected by changing several cards 
of the program of Figure 4. The response of the airframe to inputs of an 
initial velocity gust, an: initial angle of attack and a step input of elevator 
was analyzed. 

Figures 5a, 5b, and 5c show the response of the airframe velocity, 
angle of attack and pitch attitude respectively of an initial condition of 
velocity equal to ten feet per second. From this response, the phugoid 
period is found to be 33 seconds. This is in close agreement with that 
obtained in Reference 1 which was 33.5 seconds. 

Ihe airframe velocity, angle of attack and pitch attitude response to 
an initial angle of attack of two degrees is shown in Figures 6a, 6b, and 
6c. The short period mode measured from Figure 6b is 6 seconds. This 
value also agrees with that found in Reference 1 which was 6.06 seconds. 

Figures 7a, 7b, 7c, and 7d show the airframe response to a step input 
of one degree of elevator. A relationship between the angie of attack and 
elevator can be obtained from Figure 7b. A similar relationship, developed 


in Appendix III, agrees with that from Figure 7b. This relationship is: 


i = „46 & 


Aircraft With Automatic Throttle 
The response of the aircraft with the automatic throttle connected was 


investigated with the disturbance of an initial velocity gust of ten feet per 


12 


second. This disturbance was introduced with no elevator in the system and 
with the elevator function developed above as a forcing function in the system. 
The aircraft velocity and angle of attack and the engine change of thrust re- 
sponses with no elevator introduced are shown in Figures 8a, 8b, and &c. 
These same responses with the elevator input is shown in Figures 9a, 9b, 
and 9c. The use of the input of elevator as a function of angle of attack 
reduced the required change of thrust by a small amount. Use of this eleva- 
tor was made during the rest of the analysis. 

The velocity, angle of attack and change of thrust responses to an 
initial condition of two degrees to the aircraft pitch attitude is shown in 
Figures 10a, 10b, and 10c. 

The aircraft response with and without the automatic throttle as shown 
in Figures 5 through 10 agrees with the results obtained in References 1, 2, 
and 3. This indicated that the use of the digital program simulates the 
system to be analyzed as well as the analog computer utilized in the pre- 
vious investigations. 

4. Turning Flight at Constant Altitude. 

As an aircraft enters a turn and maintains an angle of bank, more lift 
must be generated in order to maintain altitude. This is due to the fact 
that the aircraft lift is no longer acting vertically, but at an angle ф , While 
the weight, of course, still acts vertically. The increased lift required can 
be generated by increasing the velocity while maintaining the same angle 
of attack or by increasing the angle of attack while maintaining the same 


velocity or by a combination of an increase in velocity and angle of attack. 
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Figure ll is a plot of the thrust required for flight for angles of bank up to 
thirty degrees. The calculations for this plot are shown in Appendix IV. 

A cross-plot of the change in thrust against angle of bank at constant 
angle of attack is plotted on Figure 12. An equation of change in thrust as 
a function of angle of bank was developed from this curve and found to be 

AT = .05ф + .855ф2 

As an aircraft maintains altitude in a tum by generating more lift, the 

normal acceleration must necessarily increase. This increase in normal 


acceleration, n, as a function of angle of bank is easily found as follows: 





을 
С08 ф = 
1 
+ عم ± اي‎ 
აპ. coso 
п = SECO -1 


In order to keep the automatic throttle from decreasing the thrust 
because of this increase in normal acceleration, the base for the normal 
acceleration loop must be changed. This can be accomplished by putting 
another input into this loop or "biasing" the loop. Figure 13 shows the 
block diagram of the system whith the expected input changes as the air- 
craft goes into a turn. The angle of attack loop can also be biased if itis 
decided to change the angle of attack in a turn. 

To investigate the use of an automatic throttle device in turning flight, 
the perturbations introduced tothethreedongitudinal equations of motion due 
to the turn must be determined. An analysis of the possible velocity and 


angle of attack perturbations is made in Appendix V. After an aircraft enters 
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a turn, if the proper perturbations are introduced and the throttle and eleva- 
tor are coordinated in such a way as to correct for these disturbances, a 
steady state should be reached at which the angle of attack, velocity and 
thrust level out at the new desired values. 

It was decided to make the turn at a constant angle of attack, thereby 
eliminating any need for biasing the angle of attack loop. This also elimi- 
nated the need to make any of the aircraft parameters variable. Utilizing 
the system shown in Figure 13, the computer program was revised to put 
the aircraft into thirty degrees angle of bank ata rate of roll of tem degrees 
per second. 

Figures 14a, 14b, and 14c show the response of the aircraft to enter- 
ing this thirty degree angle of bank turn with no perturbations introduced to 
the airframe due to the turn. When the velocity perturbation was introduced 
alone, the response obtained was that shown in Figures 15a, 15b, and 15c. 
With both the velocity and angle of attack perturbations introduced, the 
response was that shown in Figures l6a, 16b, and 16c. 

Of these three responses, the one with only the velocity perturbation 
introduced produced the results most closely agreeing with the desired 
steady state as shown in Figure 11. The response with no perturbations 
introduced shows the changes of velocity, angle of attack and thrust all 
tending toward zero after 90 seconds. When both velocity and angle of 
attack perturbations were used the response of the aircraft shows the velocity 
and angle of attack both tending toward zero while the thrust has decreased<« 


more than 9000 pounds and is still decreasing. The response of the aircraft 
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with only the velocity perturbation introduced shows the velocity steadying 
at about 19 feet per second, the angle of attack approaching zero and the 
thrust change steadying at about 1450 pounds. This is close to the desired 
steady state with the undesired being that the thrust has steadied at a value 
about 680 pounds too high. 

J: Conclusions and Recommendations. 

During this investigation the use of the digital computer was found 
very satisfactory in that originally putting the system into a useable program 
was much easier than it would have been to wire ап: analog board. Also, 
it was not necessary to do any scaling of the equations for the digital com- 
puter, as it would have been for the analog computer. И was found to be 
quite easy to make changes to the system when the digital program was used. 

The difficult part of this investigation was in simulating the aircraft 
making a turn. The perturbations introduced to the longitudinal equations 
of motion by angle of bank are not obvious and a completely satisfactory 
simulation of the turn was not realized in this investigation. A more thorough 
investigation of the effect of angle of bank on the longitudinal equations of 
motion is needed to reach a definite conclusion as to what changes must 
be made in the basic automatic throttle system to handle a turn properly. 

Another perturbation that possibly should have been included in this 
investigation is the rate of change of pitch attitude. Appendix VI develops 
a relationship between the pitch rate and the normal acceleration which is 
itself a function of the angle of bank of the aircraft. Including this pertur- 


bation may simulate the aircraft in a turn better than simulations included 
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in this report. The perturbations are introduced to the airframe equations of 


motion in the following manner. 
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It is recommended that any further investigations along these lines 
be continued as indicated in the preceding equations. It is further recom- 
mended that the possibility of using the hybrid analog-digital computer be 
investigated. This system would have the advantages of both computers and 


allow a great amount of versatility. 
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TABLE I 
F-8 Aircraft Parameters 
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Figure 1 


Block Diagram of Automatic Throttle 


Airframe System 
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Figure 3 


Flow Graph of Automatio Throttle Airframe System 
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i FIGURE 4 


COMPUTER PROGRAM OF AUTOMATIC 
THROTTLE AIRFRAME SYSTEM 


PROGRAM FLYIT 2 < 
DIMENSION Х(30)» ХООТ(30)» С(15) 
С(10) = 1.0 
PI=4e*ATANF (le) 
1 CALL INTEGI (Т»Х»ХООТ»С) 


THIS PROGRAM IS WITH AUTO-THROTTLE 
AIRCRAFT ENTERING THIRTY DEGREE ANGLE 
OF BANK TURN ы, | ს ла |- 


FOR LEVEL FLIGHT REMOVE 4 NEXT CARDS | 


® مە > 


IF (T-30) 2,393 
2 PHI=10.*T 


GO TO 4 
3 РН1-30, 
4 A=COSF(PHI*PI/180.) 
XU=-.06 
XwW=-3.3 
ZU=-.00113 
ZW=-.426 


^C ELEV*(,381/.834)*X(2) 
СЕЕ=1./А-1.+(234./32.2)*(Х(4)-Х(15)) 
XDOT(1)=XU*X(1)-32.2*X(3)+XW*X(2)+200145*X15) 
ХООТ (2) +Е|ЕУ» (- 059) +Х(4) +2И#Х(2) +2И#Х( 1) 77 
X(15)2XDOT (2) 

XDOT(3)=X(4) 

XDOT (4) =ELEV*(-2.5)-.00000455*X(5)-.338*X(4)-= 2® 
1.041%Х(15)-1.14%Х(2)%%000189%Х(1) 

ХООТ (5) + (16/1е2) # (Х(В) +Х(6)-Х(5)) 

ˆ ˆ ХООТ (6) +Х (7) ы 
XDOT(7)=-10+*X(6)-11.*X(7)-194000+*GEE 
XDOT(8)210, *C-X(8) *XC11) *XC10) &X(9)) 
X(13)=X(2) | 
XDOT(9)=-X(9)+1100.*X(13)*5703 
XDOTL10)=190.*X(13)*57.3  — 
-XDOT(11)=-2.*X(11)+1140.*X(13)%5703 

GO TO 1 : 
END 

END 


DATA CARDS MUST BE ADDED TO RUN THIS PROGRAM 
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Appendix I 


If the angle of attack 1s changed In the turn, Cy, 
Ч 
and Cp also change. The following stability derivatives 


vary with CL and Cp and ane developed below as functions 
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Appendix II 


Figure 2 18 a plot of Oy versus Cp, 
derivation'of a quadratic equation for Cp 
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Appendix III 


Using the simplified relationship between elevator 


and angle of attack given below an equation for eleva- 


tor as a function of angle can be found. 
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Appendix IV 
Below. are the formulae and calculations for the 
thrust versus velocity at various angles of attack. 


The plot of these calcutions is shown in Figure 11. 
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1. Keeping the same velocity in the turn. 
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Keeping the same angle of attack in the turn. 


“үт о 


No. 
мыл» = 241,4 Tk | see 
T= 5123 Ns. 
3. 


Flving an angle of attack corresponding to a coeffi- 


cient of lift of „8824 while tn the turn. 
me Бэх 

from Figure 2 
Cr, > . 2025 
М, = Ме де = 239,55 Tt/see 
Ка Ру, S C, 


3 


T= 585 ك‎ 


60 


&:30 
12 Keeping the same velocity in the turn. 
с 


= | 
Сыз со5% ˆ 1,003 


> m i | 
( rom F ‚aurt 7. ` р 


T = 1 pVa S Cn, | | 
T= 6666. lbs. | | 


2. Keeping the same angle of attack in the turn. 


от оа 


Үсэг = 254.5 ft/sec 





T- Y V. S CS, 
(± S560 lbs. 


3. Flying an angle of attack corresponding to a coeffi- 
cient of lift of „8824 while in the turn. 
= 8924 
Сыз. 


Съ, - ‚2025 


С, 
V NN E Е 244.5 ft/sec. 


ма ± ,ام‎ 5 С> ОЕ 


T= 50625 \bs. 


61 


Appendix V 
This appendix contatns calculations of the possible 
perturbations introduced to the velocity and the angle | 
of attack while. BÉ | 
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Utilizing the figure below, the equations following 


are east ly derived. 
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Combining the above equations the following pitch rate 


perturbat!on can be obtained. 
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